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ABSTRACT 

In  this  paper,  we  present  a method  to  estimate  the  reflected  and  fluorescent  spectral  radiance  factors  of  a fluorescent  object 
based  on  spectrophotometric  data  without  using  a monochromator.  We  use  truncated  Fourier  series  to  approximate  both  two 
spectral  radiance  factors.  Then,  based  on  the  measured  spectral  obtained  from  a spectroradiometer,  the  coefficients  of  the 
truncated  Fourier  series  are  estimated  using  an  weighted  least  square  algorithm.  The  weighting  function  is  defined  as  the  sum 
of  the  CIE  standard  x,  y,  and  z color  matching  functions.  With  the  estimated  reflected  and  fluorescent  spectral  radiance 
factors,  we  can  predict  the  color  appearance  of  a fluorescent  object  under  other  sources  such  that  the  color  difference  is 
minimized  from  viewpoint  of  human  vision. 

Keywords:  fluorescent  objects,  fluorescent  spectral  radiance  factor,  spectral  estimation,  color  appearance 

1.  INTRODUCTION 

Fluorescent  material  lias  been  widely  applied  in  various  industries.  The  fluorescent  materials  have  higher  lightness  and 
saturation  compared  with  the  traditional  non-fluorescent  materials.  The  applications  of  the  fluorescent  materials  include  (i) 
fluorescence  dyestuff,  such  as  varnish,  toys,  printery,  and  furnishing,  (ii)  fluorescent  ink,  such  as  printing,  and  warning  sign, 
and  (iii)  brightening  agents  such  as  spin,  weave,  and  papermaking.  In  industry,  the  fluorescent  materials  may  be  divided  into 
three  major  categories:  brightening  agent,  the  sunlight  fluorescent,  and  inorganic  fluorescent. 


Though  the  applications  of  the  fluorescent  materials  have  become  more  and  more  important,  it  is  difficult  to  determine 
their  appearances  under  different  sources.  Except  for  traditional  reflected  radiance,  the  fluorescent  materials  may  produce 
additional  radiance  called  fluorescence,  i.e.,  producing  low-frequency  radiant  energy  due  to  high-frequency  energy 
stimulation.  Therefore  the  chromatic  coordinates  or  the  tristimulus  values  of  fluorescent  objects  will  vary  dramatically  under 
different  light  sources.  To  correctly  predict  the  appearance  of  a fluorescent  object  under  some  specified  light  source  is  the 
main  problem  for  the  applications  of  the  fluorescent  materials.  From  the  viewpoint  of  spectrophotometry1,  the  fluorescence 
phenomenon  can  be  characterized  by  the  sum  of  the  one-dimensional  reflected  spectral  radiance  factor  /3S  and  the  two- 
dimensional  fluorescent  spectral  radiance  factor  PL  . The  key  issue  to  study  a fluorescent  material  is  to  determine  these  two 
spectral  radiance  factor  functions. 


Several  methods1,2,6  had  been  proposed  to  determine  the  spectral  reflected  radiance  factor  ps  and  spectral  luminescent 
radiance  factor  fiL  . For  the  one-monochromator  method1,  light  irradiates  on  the  fluorescent  materials  and  a monochromator 
is  used  to  control  the  wavelength  range  measured  by  detectors.  By  using  the  one-monochromator  method,  we  can  obtain  the 
overall  spectral  radiance  factor  pT  , which  is  the  sum  of  the  spectral  reflected  radiance  factor  J3S  and  spectral  luminescent 
radiance  factor  pL  . However,  the  spectral  reflected  radiance  factor  ps  and  spectral  luminescent  radiance  factor  pL  can  not 
be  distinguished  from  PT  . The  experiment  of  the  one-monochromator  method  is  easy  to  establish.  However,  as  the 
fluorescent  spectral  radiance  function  PL  can  not  be  obtained,  this  method  can  be  applied  to  determine  the  color  appearance 
of  the  fluorescent  object  under  a specified  source. 
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For  the  two-monochromator  method1 ,2,3,4,5>  two  monochromators  are  used  as  a monochromatic  radiation  selector  and  a 
radiant  detector,  respectively.  When  light  passes  through  the  first  monochromator,  the  monochromatic  radiation  at  selected 
wavelength  irradiates  the  fluorescent  material.  The  second  monochromator  then  measure  the  reflected  and  re-emitted 
radiation  from  the  fluorescent  material  at  various  wavelengths.  Therefore  we  can  obtain  a two-dimensional  reflected  matrix 
representing  the  spectrophotometric  property  of  the  fluorescent  material.  The  spectral  reflected  radiance  factor  J3S  is  usually 
determined  from  the  diagonal  of  the  matrix  and  the  spectral  luminescent  radiance  factor  pL  from  the  off-diagonal  part. 
However,  significant  error  for  the  determined  ps  may  happen  at  the  wavelengths  in  the  intersection  of  the  excitation  band 
and  the  emission  band  of  the  fluorescent  material. 

With  the  two-mode  method2,  we  can  obtain  the  overall  spectral  radiance  factor  and  spectral  conventional  reflectometer 
value  by  white  irradiating  and  monochromatic  irradiating.  Therefore,  the  overall  spectral  radiance  factor  will  equal  to  spectral 
reflected  radiance  factor  at  short-wavelengths.  And  the  conventional  reflectometer  reading  is  as  the  same  as  the  spectral 
reflected  factor  at  long-wavelengths.  Within  excitation  band  and  emission  band  of  the  fluorescent  material,  interpolation  is 
used  to  estimate  the  spectral  reflected  radiance  factor.  Hence,  the  spectral  reflected  radiance  factor  can  be  determined.  And 
the  spectral  luminescent  radiance  factor  can  be  calculated  as  the  difference  between  the  overall  spectral  radiance  factor  and 
the  spectral  reflected  radiance  factor.  Though,  we  can  determined  the  spectral  reflected  radiance,  it  is  difficult  to  obtain  a 
good  estimated  spectral  reflected  radiance  factor  in  both  excitation  band  and  emission  band  of  the  fluorescent  material. 

The  filter  reduction  method  and  the  luminescence-weakening  method  were  further  proposed2,6  to  provide  a better  method 
to  estimated  the  spectral  reflected  radiance  factor  in  fluorescent  excitation  and  emission  region.  With  the  filter  reduction 
method,  one  can  obtain  the  overall  spectral  radiance  factor  by  using  white  light  to  irradiate  the  sample.  Therefore  the  overall 
spectral  radiance  factor  will  be  equal  to  the  spectral  reflected  radiance  factor  at  short  wavelengths.  Then  the  spectral  reflected 
radiance  factor  at  long  wavelengths  will  be  obtained  by  placing  a sharp  cut-off  filter  between  the  source  and  the  sample. 
Within  excitation  band  and  emission  band  of  the  fluorescent  material,  interpolation  is  used  to  estimate  the  spectral  reflected 
radiance  factor.  A series  of  filters  are  used  to  reduce  the  fluorescent  emission  and  to  derive  an  approximation  of  the  spectral 
fluorescent  radiance  factor.  Because  the  filter  reduction  method  inducted  a series  sharp  cut-off  filter  to  give  an  approximation 
of  the  spectral  reflected  radiance  factor  within  excitation  band  and  emission  band  of  the  fluorescent  material,  this  method 
improves  estimation  accuracy  of  the  spectral  fluorescent  radiance  factor.  When  samples  have  a small  amount  of  fluorescence, 
it  is  better  to  use  the  filter-reduction  method.  On  the  contraiy,  it  is  preferable  to  use  the  luminescence-weakening  method2  in 
case  of  large  amount  of  fluorescence.  In  the  luminescence-weakening  method,  we  can  obtain  the  overall  spectral  radiance 
factor  by  using  the  white  light  to  irradiate  the  sample.  Therefore,  the  overall  spectral  radiance  factor  will  equal  to  spectral 
reflected  radiance  factor  at  short- wavelengths.  Then  the  spectral  reflected  factor  at  long-wavelengths  will  be  obtained  by 
placing  a complete  cut-off  filter  between  the  source  and  the  sample.  Within  excitation  band  and  emission  band  of  fluorescent 
material,  the  spectral  reflected  radiance  factor  will  be  deduced  by  inducting  a partly  filter  just  above  the  minimum  of  the 
overall  spectral  radiance  factor.  The  particular  wavelength  must  choose  carefully,  otherwise  it  will  induce  significant  error. 

Motta  and  Farrell7  used  sixteen  filtered  tungsten  lights  to  generate  sixteen  sources  for  which  their  spectra  can  be  used  to 
approximate  most  illuminates.  Then  a matrix  is  established  to  present  the  radiance  of  the  fluorescent  material  under  the 
sixteen  sources.  If  the  vector  representing  the  characteristics  based  on  the  sixteen  generated  sources  for  a given  source  can  be 
obtained,  then  the  colorimetric  information  of  the  fluorescent  material  under  such  a source  can  be  estimated.  However,  these 
estimated  colorimetric  data  are  not  optimized  from  human  vision. 


Based  on  the  previous  discussions,  we  now  summarize  as  follows.  For  the  the  two-monochromator  method  and  its  varieties 
such  as  the  two-mode  method,  the  filter  reduction  method,  and  the  luminescence-weakening  method  , the  reflected  spectral 
radiance  factor  ps  in  the  middle  wavelengths  is  estimated  by  interpolating  the  spectral  radiance  data  of  the  short  and  long 
wavelengths.  This  rough  interpolation  cause  estimation  error  of  ps  and  thus  of  pL  . Next,  the  monochromator  is  used  in  all 
the  work  discussed  above.  Monochromator  is  expensive  and  hard  to  maintain  and  operate.  The  most  important  thing,  which  is 
neglected,  is  that  we  should  further  validate  the  measured  or  estimated  spectral  radiance  factors,  for  example,  by  predicting 
the  color  appearance  under  some  given  source.  In  this  paper,  we  shall  construct  a method  to  estimate  the  reflected  and 
fluorescent  spectral  radiance  factor  of  a fluorescent  object  based  on  spectrophotometric  data  without  using  a monochromator. 
We  shall  use  truncated  Fourier  series  to  approximate  both  two  spectral  radiance  factors.  Then,  based  on  the  measured  spectral 
distribution,  the  coefficients  of  the  truncated  Fourier  series  will  be  estimated  using  an  weighted  least  square  algorithm.  The 
weighting  function  is  constructed  by  using  the  CIE  standard  color  matching  functions.  Then  we  shall  provide  an  estimated 
appearance  of  a fluorescent  object  under  given  sources  such  that  the  color  difference  is  minimized  from  viewpoint  of  human 
vision. 
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The  remaining  of  this  paper  is  organized  as  follows.  The  spectral  properties  of  fluorescent  materials  are  discussed  in 
Section  2.  In  Section  3,  we  shall  describe  the  setup  of  our  experimental  environment.  A least  algorithm,  which  is  used  to 
estimate  the  spectral  reflected  radiance  factor  and  the  spectral  fluorescent  radiance  factor,  and  experimental  results  are 
described  in  Section  4.  Finally,  conclusions  are  given  in  Section  5. 

2.  SPECTROPHOTOMETRIC  ANALYSIS  AND  COLORIMETRIC  MEASUREMENT  OF 

FLUORESCENT  MATERIALS 

Usually,  a fluorescent  material  is  made  of  a substrate  by  adding  some  fluorescent  agents.  Therefore,  due  to  fluorescence,  the 
overall  spectral  radiance  factor  J3r  (A,  u)  of  a material  consists  of  two  parts: 

i»rM  = AW+A(^«)  (1) 

where  fis(/ 1)  is  the  spectral  radiance  factor  reflected  by  the  substrate  and  is  that  accounts  for  the  fluorescence  effect 

of  the  material.  From  viewpoint  of  spectrophotometric  analysis  made  in  equation  (1),  to  determine  the  appearance  of  a 
fluorescent  material  irradiated  by  different  source,  it  is  important  to  find  the  two-dimensional  spectral  radiance  function 

A(A,«)  • 


Sample 


Reference  white 


Figure  1:  The  colorimetric  measurement  setup 


Now  we  turn  to  discuss  colorimetric  measurement  of  fluorescent  materials  under  the  measurement  structure  shown  in 
Figure  1.  Let  S(X) , Lprd(X),  LobJ  fiu  (A,w)  , and  Lobjsub(X)  be  the  spectral  irradiance  of  the  source,  the  spectral  radiance  of 

the  reference  white,  the  spectral  radiance  of  the  fluorescent  agents,  and  the  spectral  radiance  of  the  substrate,  respectively. 
Then  we  define  the  spectral  radiance  coefficients  of  the  reference  white  and  the  agents  of  fluorescent  sample  as 


chbj^u) : 


S(X) 


obj  _ fit 


h,U,«) 


S(A) 


(2) 

(3) 


where  u is  the  excited  wavelength  and  X is  the  emitted  wavelength.  Next,  the  spectral  radiance  factor  of  the  fluorescent 
agent  is  defined  as 


pL{X,u) 


QprdW 


(4) 


With  tire  definitions  made  in  equations  (2)  and  (3),  the  spectral  function  PL  (X,  u)  in  (4)  can  be  rewritten  as 
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(5) 


Pl(X,u) 


L0bj _ flu 


Similarly,  the  spectral  radiance  factor  of  the  substrate  of  fluorescent  material  is  defined  as 


PsW  = 


(X) 


-'obj  _ sub 


(6) 


Then  the  overall  spectral  distribution  C(A)  measured  by  the  radiometer  for  the  fluorescent  sample  is  given  by 

C(A)  = \PLQ-, u)S(u)du  + ps  (X)S(X)  (7) 

By  using  the  spectral  function  C{) l) , the  colorimetric  tristimulus  values  of  the  fluorescent  under  light  source  can  be 
calculated  as 

f C(X)x(X)dX  f C(X)y(X)ds l f C{X)~z(X)dX 

y — i no  v = i no  — 7 —ino~ 

' f S(X)y(X)dX  T f S(X)y(X)dX  ? f S(X)y(X)dX 

JX  Jx  JX 

where  x(X) , y(X) , and  z(X)  are  the  CIE  193 1 standard  color  matching  functions. 


3.  EXPERIMENT  SETUP 

3.1  Setup  of  experimental  enviroment 

In  order  to  establish  colorimetric  and  radiometric  measurements  for  the  fluorescent  samples,  we  set  up  an  experiment  as 
shown  in  Figure  1.  The  measurement  instrument  setup  follows  the  45/0  geometry.  Both  the  fluorescent  sample  and  the 
reference  white  are  irradiated  by  a source  and  their  spectral  radiance  are  measured  by  a spectroradiometer. 

3.2  Selection  of  Light  Sources 

To  select  the  sources,  we  must  consider  the  fluorescence  phenomena  of  the  samples.  First,  the  fluorescence  is  excited  within 
the  short-wavelength  range  of  the  incident  radiant  energy.  Next,  the  intensity  of  the  fluorescent  spectral  radiance  factor 
PL  (A,  u)  is  less  than  that  of  the  reflected  spectral  radiance  factor  PS(X) . Therefore  a source  with  its  power  distributed  over 
long-wavelength  range  can  be  used  to  effectively  estimate  the  spectral  function  ps  ( X ) . However,  to  measure  the  fluorescent 
spectral  function  PL  (X,  u) , we  need  to  select  a source  with  rich  excitation  over  short-wavelength  range.  We  use  seven 
different  sources  in  our  experiment.  These  sources  are: 

Source  1:  a Tungsten  lamp, 

Source  2:  a test  lamp  from  CVI  corporation, 

Source  3:  a Xenon  lamp, 

Source  4:  an  approximate  of  the  illuminant  D65  , 

Source  5:  the  first  test  fluorescent  lamp  , 

Source  6:  the  second  test  fluorescent  lamp. 

Source  7:  the  third  test  fluorescent  lamp. 

The  relative  spectral  power  distribution  are  show  in  Figure  2,  Figure  3,  Figure  4,  Figure  5,  Figure  6 , Figure  7,  and  Figure 
8,  respectively.  Note  that  the  power  spectrums  of  Source  1 and  Source  2 are  mainly  distributed  over  long-wavelength  range 
while  the  excitation  of  short-wavelength  range  is  richer  for  Source  3,  Source  4,  Source  5,  Source  6,  and  Source  7. 

3.3  Descriptions  of  the  Fluorescent  Samples 

We  will  estimate  the  color  appearance  of  three  kinds  of  fluorescent  samples  under  different  sources.  There  are  three 
fluorescent  samples,  which  are  yellow,  pink,  and  white,  respectively,  to  be  tested. 
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Releative  spectral  power 


Source  1 Source  3 


Figure  2:  Relative  spectral  power  distribution  of  Source  1 
with  2 run  sampling  between  380-780nm 


Source  2 


Figure  3 : Relative  spectral  power  distribution  of  Source  2 
with  2nm  sampling  between  380-780nm 


Figure  4:  Relative  spectral  power  distribution  of  Source  3 
with  2nm  sampling  between  380-780nm 


Source  4 


Figure  5:  Relative  spectral  power  distribution  of  Source  4 
with  2nm  sampling  between  380-780nm 
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Source  5 


Source  7 


Figure  6:  Relative  spectral  power  distribution  of  Source  5 
with  2nrn  sampling  between  380-780nm 

Source  6 


Figure  8:  Relative  spectral  power  distribution  of  Source  7 
with  2mn  sampling  between  380-780nm 


Figure  7:  Relative  spectral  power  distribution  of  Source  6 
with  2nm  sampling  between  380-780nm 

3.4  The  Observer 

The  definition  of  observer  in  our  experiment  follows  the  CIE  1931  Standard  Observer.  The  purpose  of  using  the  45/0 
geometry  is  to  avoid  re-excitation  of  the  fluorescent  materials.  We  use  a spectroradiometer  (PR-704)  to  measure  the  spectral 
distribution  of  the  sample.  The  spectral  distribution  is  sampled  every  2nm  within  the  wavelength  range  from  380nm  to 
780nm. 

3.5  Measurement  Procedure 

The  spectral  radiance  of  the  reference  white  irradiated  by  the  three  different  sources,  Source  1 to  Source  3,  are  measured  by 
the  spectroradiometer.  Then,  the  spectral  radiance  function  C(A)  for  each  tested  fluorescent  sample  is  measured  with  respect 
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to  the  three  light  sources.  In  order  to  reduce  the  effect  of  random  fluctuation1,  each  spectral  distribution  is  obtained  by 
averaging  ten  measurements. 


4.  SPECTRAL  ESTIMATION  ALGORITHM  AND  EXPERIMENTAL  RESULTS 

To  estimate  the  spectral  functions  ps (A)  and  /?L(A,w) , we  use  Fourier  series  expansions  to  represent  these  two  functions, 
and  equation  (7)  is  transformed  into  the  linear  regression  form  in  order  to  apply  standard  least  square  algorithm.  In  particular, 
to  minimize  the  estimation  error  of  the  color  appearance  of  a fluorescent  sample,  the  cost  function  involved  in  the  least 
square  algorithm  is  weighted  by  a function  which  is  the  summation  of  the  CEE  color-matching  functions  3c(A) , y(A)  , and 

z(A) . For  each  fluorescent  sample,  the  measured  spectral  distribution  under  Source  1,  Source  2,  and  Source  3,  are  used  in  the 

least  squared  algorithm  to  obtain  the  estimated  functions  ps  (A)  and  PL  (A,  u) . These  estimated  spectral  radiance  functions 
are  then  used  to  predict  the  tristimulus  of  the  tested  fluorescent  samples  under  Source  4,  Source  5,  Source  6,  and  Source  7. 


4.1  Least  Square  Estimation  Algorithm 


For  representing  the  spectral  function  PL  (A,  u) , we  use  a truncated  two-dimension  Fourier  sine  series  expansion;  and  for  the 
spectral  reflected  radiance  factor  ps  (A) , a truncated  one-dimension  Fourier  sine  series  expansion  is  used  as  follows.  The 
estimated  fluorescent  spectral  radiance  factor  is  of  the  form: 


- . x • ,w(A-A0)fl\  . w(w-w0);r 

Pi  (A,  «)  = XI  6m  sin(  0 -)  sin(  — ; ■—■) 


K 


(8) 


where  A=Z,=400nm,  A0  = uQ  = 380 nm , and  Qnm  are  the  coefficients  to  be  determined.  The  constant  N-26  is  the  number  of 
sinusoidal  function  used  in  the  Fourier  series  expansion  with  respect  to  both  the  u axis  and  the  A axis.  The  structure  of  the 
estimated  reflected  spectral  radiance  function  ps  (A)  is  given  by 


Aa)=IX  sin(^3^W0 

n=\ 


(9) 


where  =400nm  and  S =51  is  the  number  of  coefficients,  8a  and  8Q , in  this  series  to  be  determined.  Note  that  the  number 
of  all  the  coefficients  to  be  determined  in  (8)  and  (9)  is  676.  This  number  is  much  less  than  that  required  in  the  two- 
monochomator  method. 


Since  the  spectral  distribution  C( A)  is  measured  every  2nm,  the  estimated  total  spectral  radiance  distribution  C(A)  is 
defined  as 


- _ S vfr!.  /i(A-A0);r  . m(u-u0)n  ^ 

C(A)  = XZXms,n(— — --  -)  sm( — — ; — — )S(«)A«+ 


w-380  n=\  m=l 
u even 


K 


IX  sin (g^oj£)+<?o  m)  (10) 


where  Au  = 2 nm  . The  above  expression  can  be  rewritten  into  linear  regression  which  fits  the  setup  in  least  square 
estimation.  Define  a regression  vector  h(A)  and  a parameter  vector  <f>  as 

hW  = U.i  W ■ ■ ■ /i,w-i  (A)  fi,\  (A)  • • • /a,.,,*.,  (A)  r,  (A)  • ■ ■ r9_,  {A)  r0  (A)f 
$ - 1^1,1  ®1,\  '"@n-\,n-\  ^oY 

where  the  superscript  T denotes  matrix  transpose,  r0  (A)  = S(A)  , and 

( \ 


, . ,n(A-A0)n 

f nm  (A)  = sin( 

K 


780 

z 

«=380 


vw 


n,  m = 1, . . . , N - 1 


ra  (A)  = sin( 


. a(A-Aa)x 


O 


)S(A),  a = \,...S-\ 
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(11) 


Then  equation  (10)  can  be  rewritten  into  a linear  regression  form  as 

C(X)=hT  {A)<f> 

For  any  given  estimated  vector  <p  , the  spectral  estimated  error  at  wavelength  A is  given  by  C(X)-hT  . Let  £(X)  be  a 
weighting  function  defined  as 

§(X)  = x(X)+y(X)  + z(X)  (12) 

where  x{X)  , y(X) , and  z(X)  are  the  CEE  1931  standard  color-matching  functions.  The  weighting  function  g(A)  is  used  to 
emphasize  the  effective  wavelength  range  for  human  vision.  By  the  definitions  of  the  CEE  1931  standard  color-matching 
functions,  we  can  find  that  the  effective  wavelength  region  for  human  vision  is  between  380nm  and  750nm  which  is  the 
union  of  the  supports  of  the  functions  x(X) , y(X)  „ and  z(X) . In  order  to  minimize  the  color  difference  from  human  vision 

coordinate  (X,  Y,  Z),  the  weighting  function  £(2)  is  included  in  the  least  square  algorithm  to  specify  the  wavelength  range 
within  which  the  spectral  estimation  error  is  more  significant. 


There  are  three  sets  of  spectral  distributions  irradiated  by  Source  1 to  Source  3 to  be  used  for  training  the  estimated 
parameter  vector  (j) . Each  spectral  set  is  sampled  every7  2nm  from  380nm  to  780nm.  Denote  Cj  (X),  j = 1,2,3,  as  the 
measured  overall  spectral  radiance  function  and  hj(X)  as  its  related  regression  vector  under  Source  1,  Source  2,  and  Source 
3,  respectively.  Define  a sequence  of  wavelengths  as  Xj  =380  + 2 j ran,  j = 0, . . . ,200  . Also  define  the  following  vectors  and 
matrix  as 

V = to-^oo)  ^a0)-^200)  #(A0)-«<lMo)] 

C_T  = [C,(A0)-C,(A200)  C2(l0)-C2(l200)  C3(A0)-C3(A200>] 
hT  = [//,(A0  )•••*!  (^joo ) h2(A0)-h2(Am)  h3(A0)---h3(A200  )] 

Now  denote  £.  and  C.  as  the  i-th  entry  of  the  vectors  and  C , respectively.  Let  h{  be  the  transpose  of  the  z'-th  row  vector 
of  h . Then  the  cost  function  to  be  minimized  in  the  least  square  algorithm  is  defined  as 

1 603  _ 2 

(13) 

2 M 

The  best  parameter  vector  estimate  which  minimizes  the  cost  function  J{(p)  is  denoted  as  (j)LS  , i.e. , 

< Ks  = arg  min  J (<j)) 


A recursive  algorithm  to  compute  <j)LS  is  given  by 


<t>t  = 


p^  = 


1 +hT,-iPi_2hi_l 


l + /7r,-lP._2/zM 


i - 2, ...,603 


(14) 


where  the  initial  data  <fix  can  be  arbitrarily  given  and  P0  must  be  a positive  definite  matrix.  The  last  estimate  <j)m  is  set  to 

$LS  • 

By  the  proposed  weighted  least  square  algorithm,  we  can  obtain  the  estimated  reflected  spectral  radiance  function  ps  (X) 
and  the  estimated  fluorescent  spectral  radiance  function  PL(X,u)  for  each  fluorescent  sample.  By  using  equation  (10),  we 

can  construct  the  spectral  radiance  function  C(X)  that  estimates  the  spectral  radiance  function  C(X)  for  the  test  fluorescent 
samples  under  Source  4,  Source  5,  Source  6,  and  Source  7.  Comparisons  of  the  estimated  overall  spectral  radiance  function 
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C(/ 1)  and  the  measured  one  C(A)  for  the  three  fluorescent  samples  under  different  sources  are  illustrated  in  Figure  9 to 
Figure  14.  In  these  figures,  the  measured  spectral  function  C(A)  is  well  approximated  by  its  estimated  one.  The 
approximation  error  for  the  measured  spectral  radiance  function  CQ i)  leads  to  color  difference.  Color  differences  in  terms  of 
CIELAB  Color  Difference  Formulae  for  each  fluorescent  sample  under  Source  1 to  Source  7 are  presented  in  Table  1.  The 
values  of  A E*b  under  Source  1,  Source  2,  and  Source  3 in  Table  1 are  very  small  because  the  measured  spectral  distribution 
under  these  sources  are  used  to  train  the  parameter  estimates  in  the  least  square  algorithm.  The  estimated  spectral  functions 
(A)  and  fiL  ( A , u)  are  then  used  to  predict  the  color  appearance  of  fluorescent  samples  under  Source  4,  Source  5,  Source 
6,  and  Source7.  Since  the  measured  spectral  data  under  Source  4,  Source  5,  Source  6,and  Source7  are  not  used  in  the  least 
square  algorithm,  the  values  of  A E*ab  under  these  sources  are  usually  larger.  It  is  more  difficult  to  decrease  the  values  of 

A E*ah  under  Source  5,  Source  6,  and  Source  7 since  these  sources  are  mainly  based  on  fluorescent  lamps,  which  have  discrete 
spectral  distributions. 


Source  1 

Source  2 

Source  3 

Source  4 

Source  5-7 

Sample  1 

0.0010 

0.0001 

0.0011 

3.55 

6. 73 (under  Source  5) 

Sample  2 

0.0013 

0.00001 

0.0014 

i 

3.59 

8. 81  (under  Source  6) 

Sample  3 

0.0130 

0.0060 

0.0003 

4.6 

4.41  (under  Source  7) 

Table  1 : Color  differences  of  the  fluorescent  samples  under  different  sources 


5.  CONCLUSIONS 

In  this  paper,  we  have  presented  a method  to  estimate  the  reflected  and  fluorescent  spectral  radiance  factor  of  a fluorescent 
object  based  on  spectrophotometric  data  without  using  a monochromator.  We  use  truncated  Fourier  series  to  approximate 
both  two  spectral  radiance  factors.  The  total  number  of  parameters  to  be  determined  and  recorded  are  much  lees  than  that  in 
the  two-monochromator  method.  With  the  special  weighted  least  square  algorithm,  we  are  able  to  provide  estimated  color 
appearance  of  test  fluorescent  objects  under  different  sources  such  that  the  color  difference  is  minimized  from  viewpoint  of 
human  vision.  With  the  experiment  test  made  in  Section  4,  the  proposed  method  can  be  effectively  used  to  predict  color 
appearance  for  tire  applications  of  fluorescent  materials  such  as  fluorescent  dyestuff,  fluorescent  ink,  and  brightening  agents. 
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Figure  9:  Comparison  of  C(X)  and  C(X)  for  the  yellow 
fluorescent  sample  under  Source  4. 
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Figure  10:  Comparison  of  C(A)  and  C(X)  for  the  pink 
fluorescent  sample  under  Source  4. 
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